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ABSTRACT 

Using high signal-to-noise echelle spectra of 3C 273 obtained with the Space Tele- 
scope Imaging Spectrograph (resolution = 7 km s" 1 FWHM), we constrain the metal- 
licities of two Lya clouds in the vicinity of the Virgo cluster. We detect C II, Si II, and 
Si III absorption lines in the Lya absorber at z^ s = 0.00530. Previous observations 
with the Far Ultraviolet Spectroscopic Explorer have revealed Ly/3— Ly6* absorption lines 
at the same redshift, thereby accurately constraining the H I column density. We model 
the ionization of the gas and derive [C/H] = -1.2+ ,; 3 ,, [Si/C] = +0.2 ± 0.1, and log 
71-h = —2.8 ± 0.3. The model implies a small absorber thickness, ~ 70 pc, and thermal 
pressure p/k « 40 cm -3 K. It is most likely that the absorber is pressure confined by an 
external medium because gravitational confinement would require a very high ratio of 
dark matter to baryonic matter. Based on a sample of Milky Way sight lines in which 
carbon and silicon abundances have been reliably measured in the same interstellar 
cloud (including new measurements presented herein), we argue that it is unlikely that 
the overabundance of Si relative to C is due to depletion onto dust grains. Instead, this 
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probably indicates that the gas has been predominately enriched by ejecta from Type II 
supernovae. Such enrichment is most plausibly provided by an unbound galactic wind, 
given the absence of known galaxies within a projected distance of 100 kpc and the 
presence of galaxies capable of driving a wind at larger distances (e.g., H I 1225+01). 
Such processes have been invoked to explain the observed abundances in the hot, X-ray 
emitting gas in Virgo. However, the sight line to 3C 273 is more than 10° away from 
the X-ray emission region. We also constrain the metallicity and physical conditions of 
the Virgo absorber at z a b s = 0.00337 in the spectrum of 3C 273 based on detections of 
O VI and H I and an upper limit on C IV. If this absorber is collisionally ionized, the 
O Vl/C IV limit requires T > 10 53 K in the O Vl-bearing gas. For either collisional 
ionization or photoionization, we find that [O/H] > —2.0 at z a b s = 0.00337. 

Subject headings: galaxies: clusters: individual (Virgo) — ISM: abundances — inter- 
galactic medium — quasars: absorption lines — quasars: individual (3C 273) 

1. Introduction 

The heavy element enrichment of the intergalactic medium (IGM) provides a valuable con- 
straint on physical processes that play critical roles in the evolution of galaxies. Elements heavier 
than boron are synthesized in stars (presumably mostly located within galaxies), and understanding 
how these elements are removed from galaxies and transported into the IGM will likely yield insights 
on galaxy evolution. Is the IGM enriched by supernova-driven galactic winds, dynamical processes 
such as tidal or ram-pressure stripping, or all of these mechanisms? How does this affect /regulate 
star formation in galaxies? IGM enrichment also has cosmological implications. For example, the 
processes which deposit heavy elements in the IGM will also deposit energy, and this "feedback" 
can have important effects on cosmological structure growth (e.g., Kaiser 1991; Ponman, Cannon, 
& Navarro 1999; Cen & Bryan 2001; Voit & Bryan 2001). 

IGM enrichment has been modeled using methods ranging from semi-analytic calculations to 
hydrodynamic simulations, and a variety of theoretical predictions have been made regarding the 
heavy element content of intergalactic gas at both high and low redshifts (e.g., Cen & Ostriker 
1999b; Ferrara, Pettini, & Shchekinov 2000; Aguirre et al. 2001). These predictions can be tested 
using QSO absorption lines to estimate the metallicity of the IGM in a variety of contexts. Observa- 
tions of QSO absorption-line systems have shown that the high-redshift IGM is remarkably-widely 
enriched with metals (Cowie et al. 1995; Tytler et al. 1995), and statistical analyses have indicated 
that this enrichment extends even to the absorbers with the lowest H I column densities (Ellison 
et al. 2000; Schaye et al. 2000). 

Recently, high signal-to- noise (S/N) UV spectroscopy of 3C 273 with the Space Telescope 
Imaging Spectrograph (STIS) and the Far Ultraviolet Spectroscopic Explorer (FUSE) has provided 
an opportunity to probe the IGM enrichment within a particularly interesting environment: the 
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Virgo cluster/supercluster. In this paper we derive constraints on the heavy element abundances in 
two Lya clouds in the immediate vicinity of this cluster. The sight line to 3C 273 extends through 
the Virgo cluster "Southern Extension" (e.g., Binggeli, Popescu, &; Tammann 1993), although the 
projected distance from the sight line to M87 and the X-ray region is substantial (> 3 Mpc). In the 
3C 273 spectrum, there are three Lya clouds at Virgo redshifts; in this paper we are mainly inter- 
ested in the Virgo absorbers at z a bs = 0.00337 and 0.00530. The connections between these Virgo 
absorption systems and nearby galaxies have been extensively studied (see §4.3), but measurements 
of the metallicities of these Lya absorbers have not been previously obtained. It is now possible 
to estimate the metal abundances in these Virgo Lya clouds due to two observational advances. 
First, Sembach et al. (2001) have presented a surprising result from the FUSE observation: the 
Virgo absorber at z a b s = 0.00530 is detected in seven higher Lyman series lines (Ly/3 - Ly0). This 
clearly establishes that the H I column is substantially higher than previously thought; from a 
curve-of-growth (COG) analysis, Sembach et al. (2001) derive log N(H I) = 15.85toQg, an increase 
by a factor of more than 40 over the previous estimate based on fitting the profile of Lya only. 11 
More importantly, the H I column density and Doppler parameter are accurately constrained from 
the COG. Second, the high signal-to-noise STIS observations at 7 km s _1 resolution have enabled 
the detection of very weak lines in the 3C 273 spectrum, including weak lines of C II, Si II, and 
Si III at z a bs = 0.00530. In this paper, we use these weak metal lines and the new H I measure- 
ment to explore the heavy element enrichment of this Virgo Lya system at z a b s = 0.00530. We 
also constrain the metallicity of the lower-redshift Virgo Lya cloud (z a b s = 0.00337) based on the 
detection of O VI and the absence of C IV. 

The paper is organized as follows. In §2 we present the new STIS observations and absorption 
line measurements, and we briefly comment on the well-known galaxy structures at Virgo redshifts 
in the vicinity of the 3C 273 sight line. In §3 we model the ionization of the gas and derive constraints 
on the metallicity, density, and temperature of the Lya absorbers. We discuss the implications of the 
measurements and the nature of the absorbing gas in §4, and we summarize the main results in §5. 
In the Appendix we provide additional details on new Galactic ISM measurements to supplement 
the analysis in §4. Throughout the paper, we report heliocentric wavelengths and redshifts, 12 and 
we assume Hq = 75 km s _1 Mpc -1 and go = 0. 



11 Using ORFEUS //observations of the Ly/3 line at this redshift, Hurwitz et al. (1998) provided the first evidence 
that the column density of this system was substantially higher than expected based on profile fitting of the Lya 
transition alone. However, with only Lya and Ly/3, they were unable to tightly constrain 7V(H i) and 6 (see their 
Figure 4). 

12 In the direction of 3C 273, i>lsr = vq + 2.3 km s _1 assuming the standard definition of the Local Standard of 
Rest (Kerr & Lynden- Bell 1986). 
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2. Observations and Measurements 

Observations of 3C 273 were obtained with the E140M echelle mode of STIS, using the C/2 x $'.2 
aperture, on 2000 May 2 and 2000 June 21-22. This mode provides a resolution of 7 km s" 1 (FWHM) 
and records the spectrum between 1150 and 1700 A with only a few small gaps between orders at 
A > 1615 A (Kimble et al. 1998; Woodgate et al. 1998). The total integration time was 18.67 ksec, 
and this resulted in signal-to-noise ratios ranging from ~6 to 66 per resolution element. The data 
reduction followed standard procedures, including the scattered light correction developed by the 
STIS Investigation Definition Team. 

A selected portion of the final spectrum is shown in Figure 1. This portion of the spectrum 
spans the wavelength range of Lya lines at redshifts within the Virgo cluster/supercluster. Previous 
HST studies have detected two Lya clouds within the Virgo cluster at high significance levels at 
z a bs = 0.003 and 0.005, but a third Lya absorber at = 0.007 was reported in some studies but 
not confirmed (see §3.3.1 of Morris et al. 1993). For example, Brandt et al. (1997) do not list a 
line at z a b s = 0.007 in their intermediate-resolution GHRS atlas of 3C 273. However, Penton et 
al. (2000) claim a significant detection of this line using GHRS/G160M data, and a weak feature 
is evident at the right wavelength in the Brandt et al. spectrum. The new STIS spectrum settles 
this question: a broad feature with rest-frame equivalent width 13 W T = 94 ± 14 mA is detected at 
z abs = 0.00745. 

Figure 1 also indicates the expected wavelengths of Lya lines at the mean redshifts of several 
Virgo galaxy structures near the 3C 273 sight line. It has long been recognized that there are 
several large "clouds" of galaxies in the immediate vicinity of the Virgo cluster that have distinct 
kinematics and morphologies (e.g., de Vaucouleurs 1961). According to Binggeli, Popescu, & 
Tammann (1993), the 3C 273 sight line pierces the structures known as the Virgo "Southern 
Extension" (< v >= 1012 ± 51 km s _1 , velocity dispersion a = 198 km s" 1 ) and the Virgo W 
cloud (< v >= 2099 ± 41 km s -1 , a = 356 km s _1 ). The Virgo W cloud is also in this region 
(< v >= 1310 ± 80 km s -1 , a = 301 km s _1 ). The relationships between these structures and the 
Virgo cluster proper are not entirely clear. As noted, our pencil beam to 3C 273 is at a substantial 
projected distance from the Virgo X-ray emission. At any rate, the sight line probes a structure that 
contains a large number of galaxies in the general vicinity of the Virgo cluster. For convenience, 
we simply refer to this region as the Virgo supercluster in this paper. 

Figure 2 shows the absorption profiles of lines detected in the z a b s = 0.00530 absorber including 
H I Lya, Si II A1260.4, Si III A1206.5, and C II A1334.5, as well as the regions of the spectrum 
covering the undetected Si IV A1393.8 and C IV A1548.2 lines. We note the presence of a very 
weak feature near the expected wavelength of the Si IV transition. However, the significance of this 
feature is much less than 3a, and furthermore this spectral region is partially obscured by the field 
electrode of the FUV-MAMA detector, the "repeller wire" (Woodgate et al. 1998). In the majority 



Rest-frame equivalent width W T = Wobs/(l + «abs), where Wobs is the observed equivalent width. 
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Fig. 1. — Portion of the STIS echelle spectrum of 3C 273 plotted versus observed heliocentric 
wavelength. In this figure, the spectrum has been binned (two pixels into one) for display purposes 
only (all other figures and measurements make use of the full-resolution, unbinned data). The three 
Lya clouds within the Virgo cluster region at z a b s = 0.00337, 0.00530, and 0.00745 are indicated. 
The redshifted Lya wavelengths corresponding to the mean velocities of the "Southern Extension" , 
W, and W galaxy clouds in Virgo (from Binggeli et al. 1993) are also indicated; the horizontal 
bars indicate the velocity dispersions of these Virgo structures. The solid line near zero is the la 
flux uncertainty, adjusted to reflect the binning. The broad absorption feature centered at 1216 A 
is the damped Lya line due to H I in the Milky Way ISM, and the spike is the geocoronal Lya 
emission line. 
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Fig. 2. — Absorption profiles of the H I Lya, Si II A1260.4, Si III A1206.5, and C II A1334.5 lines associated with 
the Lya absorber at z a hs = 0.00530 in the 3C 273 spectrum, plotted vs. rest-frame velocity where v = km s _1 at 
Zabs = 0.00530. The spectral regions of the Si IV A1393.8 and C IV A1548.2 lines, which are not detected, are also 
shown. The dotted lines show the continua adopted for the measurements and upper limits listed in Table 1. Note 
that the expected wavelength of the Si IV line is near the region affected by the FUV MAMA repeller wire, and the 
feature at v ~ 300 km s _1 is an artifact due to the repeller which was not adequately removed by the flat field. The 
smooth solid curve in the top panel shows a theoretical Lya profile with log AT(H I) = 15.85 and b — 16.1 km s _1 , as 
derived by Sembach et al. (2001) from a curve-of-growth analysis of FUSE observations of 3C 273. The line spread 
function from the STIS Handbook has been assumed to produce this theoretical profile, and its velocity centroid has 
been adjusted to best fit the observed line. The observed Lya profile is ~ 50 mA stronger than expected; the excess 
absorption is evident primarily in the blue wing of the line. 
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Fig. 3. — Absorption profile of the H I Lya line at z a b s — 0.00337 in the direction of 3C 273, plotted 
versus rest-frame velocity where v = at z a ^ s = 0.00337. The spectral regions of the undetected 
Si III A1206.5 and C IV A1548.2 lines are also shown. 
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of the 3C 273 observations, the repeller introduces significant artifacts into the spectrum which are 
not adequately removed by the flat field. Sembach et al. (2001) do not report significant detection 
of any metals at z a b s = 0.00530 in the FUSE 3C 273 spectrum, but the stronger transition of the 
O VI doublet is blended with strong Milky Way absorption lines. Similarly, the strong C III A977.0 
line falls close to a Galactic H 2 absorption line at this redshift. 

The only line detected at z abs = 0.00337 in the STIS E140M spectrum of 3C 273 is the Lya 
transition. However, Sembach et al. (2001) have also detected Ly/3 and O VI A1031.9 at this 
redshift in the FUSE spectrum. Figure 3 shows the Lya line at z a b s = 0.00337 as well as regions 
of strong metal lines that are not detected. 

We have measured the absorption-line equivalent widths and integrated apparent column densi- 
ties 14 using the techniques of Sembach & Savage (1992), which include in the overall error evaluation 
the uncertainty in the level and curvature of the continuum and a 2% flux zero point uncertainty, as 
well as the usual photon-counting uncertainty. The equivalent widths and column densities of the 
Zabs = 0.00530 absorption system are summarized in Table 1 along with upper limits on undetected 
species of interest. Equivalent widths, column densities, and upper limits from the z a bs — 0.00337 
absorber are listed in Table 2. 

Based on previous GHRS observations of 3C 273, Sembach et al. (2000) report that the 
equivalent width of the Lya line at z a b s = 0.00530 is ~40 mA stronger than predicted by their 
COG analysis of the Ly/3 - Ly6 lines at this redshift. As shown in the top panel of Figure 2, we 
confirm this result with the higher resolution STIS data. The excess equivalent width in the STIS 
spectrum is ~50 mA and is most evident on the blue side of the profile. While this establishes that 
multiple components contribute to the H I absorption, the additional component in the Lya profile 
is probably too weak to significantly affect the lines higher in the series, and Sembach et al. show 
that including or excluding Lya does not change N(R I) and b significantly (see their Table 4). We 
discuss the implications of the Lya profiles below. 



14 The integrated apparent column density is given by iV a = (m c c/ire 2 )(f\)~ 1 f T a (v)dv — 3.768 x 
10 14 (/A)~ 1 J \n[I c (v)/I(v)]dv, where / is the oscillator strength, A is the wavelength (the numerical coefficient in 
the latter expression assumes A is in A), I(v) is the observed line intensity and I c (v) is the estimated continuum 
intensity at velocity v; see Savage & Sembach (1991) for further details. 
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Table 1. Equivalent Widths and Integrated Column Densities of the Virgo Lya Absorber at z a b s 

= 0.00530 a 



Species 


Ao b 
(A) 




(mA) 


< v > d 
(km s _1 ) 


log N 


S/N e 


HI 


1215.67 


0.416 


377±12 


-10 ±3 


> 14.3 


23 


Sill 


1260.42 


1.18 


8.9±2.6 


-12 ±5 


11 76+ 011 

LL - ' D -0.15 


30 


Si III... . 


1206.50 


1.67 


36.8±6.1 


-5 ±2 


12.33±0.08 


11 - 16 


SiIV.... f 


1393.76 


0.514 


< 8.7§ 




< 12. h 


27 


CII 


1334.53 


0.127 


10.3±2.6 


-9 ±3 


12.74;°;}° 


34 


CIV 


1548.20 


0.191 


< 12.0 g 




< 12.5 h 


25 


Al II 


1670.79 


1.83 


< 17.3 g 




< 11.6 h 


17 


Fe II 


1608.45 


0.0580 


< 25.3 




< 13.3 


16 


I 


1302.17 


0.0519 


< 6.7 s 




< 12.9 h 


34 



a With the exception of Lya, all quantities are integrated from —30 to 20 km s 1 where 
v = km s _1 at z a b s = 0.00530. The Lya quantities are integrated from —160 to 70 km 

S ~\ 

b Rest-frame vacuum wavelength and oscillator strength from Morton (2002) or Morton 
(1991). 

c Rest-frame equivalent width. 

d Profile- weighted mean velocity of the line, < v >= J v[l — I (v) / I c (v)]dv / J[l — 
I(v)/I c (v)]dv. 

c Signal-to-noise ratio per two-pixel resolution element in the continuum near the absorp- 
tion line. The Si III A1206.5 line is located in the blue wing of the Milky Way damped 
Lya profile. Here the steep slope of the continuum leads to rapidly decreasing S/N with 
increasing wavelength, so we list the S/N range within ±50 km s _1 of the line. 

f The Si IV A1393.8 line is located in a region of the spectrum affected by the FUV 
MAMA repeller wire. This introduces features into the spectrum which are not adequately 
removed by the flatfield. Furthermore, the line is close to the order edge where spurious 
features are occasionally evident. 

s 3er upper limit. 

h 3cr upper limit derived from the equivalent width upper limit assuming the linear curve- 
of-growth is appropriate. 
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Table 2. Equivalent Widths and Integrated Column Densities of the Virgo Lya Absorber at z a b s 

= 0.00337 a 

Species A / W r <v> log N S/N 

(A) (mA) (km s _1 ) 

HI 1215.67 0.416 389±12 -3 ± 1 > 14.2 17-20 

Si III 1206.50 1.67 <34.5 b ••• <12.2 b 20 

CIV 1548.20 0.191 < 27.6 b ••• < 12. 8 b 24 

OVI 1031.93 0.0659 25.4±7.1 c ••• 13.32l°;^ c 

a All quantities are integrated from —100 to 100 km s _1 where v = 
km s _1 at z a bs = 0.00337, except the O VI measurements, which are from 
Sembach et al. (2001). See Table 1 footnotes for further details. 

b 3o" upper limit. 

c Equivalent width and column density reported by Sembach et al. (2001). 
The equivalent width is the weighted mean of the measurements from the 
two FUSE channels, and the column density is from their Table 4. 
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3. Analysis 
3.1. z abs = 0.00530 

From the COG analysis of the Ly/? - Ly# lines at z a b s = 0.00530, Sembach et al. (2001) obtain 
b = 16.1 ± 1.1 km s -1 . Therefore the gas temperature T < 16,000 K. Even though iV(H I) has 
been dramatically revised to a larger value in this absorber, with log iV(H I) = 15.85 the gas is still 
optically thin below the Lyman limit. With this column density, the gas is likely to be substantially 
ionized by the UV background radiation from QSOs. Consequently, to derive abundances, ionization 
corrections must be applied. For this purpose, we have made use of ionization models constructed 
with the photoionization code CLOUDY (v94.00, Ferland et al. 1998). Collisional ionization 
may also be important, but collisional processes are included in CLOUDY, and as we shall see, 
the mean gas temperature in the best CLOUDY model is fully consistent with the upper limit 
implied by the b— value. Therefore even if collisional ionization is dominant, these models should 
provide reasonable estimates of the metal abundances. Throughout this section we shall restrict 
our attention to gas-phase abundances; we discuss how depletion onto dust grains could affect the 
abundances in §4. 

In the photoionization models, the gas is assumed to be adequately approximated by a plane- 
parallel, constant density slab. We adopt the background radiation field at z = due to QSOs, as 
computed by Haardt & Madau (1996) with the mean intensity at 1 Rydberg set to J u = 1 x 10~ 23 
ergs s _1 cm~ 2 Hz -1 sr _1 (see Shull et al. 1999; Dave & Tripp 2001; Weymann et al. 2001 and 
references therein). 15 The model column densities primarily depend on the ionization parameter 
U or, equivalently, the gas density since U = n^/nu = H ionizing photon density/total H number 
density, and n 7 is set by the assumed shape and intensity of the ionizing radiation field. At low 
metallicities, the column densities scale directly with the overall metallicity Z. We constrain the 
ionization parameter using the observed Si Ill/Si II ratio, but we allow the relative abundance of 
silicon to carbon to vary since the gas may not have a solar abundance pattern. 

We begin with the most simple model, assuming that the C II, Si II, and Si III absorption lines 
arise in the same gas. As we discuss below, there is some evidence that these lines do not all occur 
in the same phase, but the evidence is not yet strong enough to rule out this possibility. With 
this assumption, Figure 4 shows the CLOUDY model that provides the best simultaneous fit to 
the observed C II, Si II, and Si III column densities. In the usual logarithmic notation, 16 we obtain 
an excellent fit with the following abundances: [C/H] = —1.4 and [Si/H] = —1.0 with «h ~ 10 -31 
cm~ 3 (or log U ~ —3.35). At this density, the model predicts a gas temperature of ~ 14,000 K, 



15 Scott ct al. (2002) have recently derived J v = 6.5±^ 8 6 x 10~ 23 ergs s _1 cm~ 2 Hz -1 sr _1 from the proximity effect 
in a large number of QSO spectra obtained with the Faint Object Spectrograph at 0.03 < z < 1.0. Given the redshift 
range to which this applies, this is consistent with the expected evolution of J„ and the value we have assumed for 

2^0. 

16 [X/Y] = log (X/Y) — log (X/Y) Q . Solar reference abundances are taken from Holweger (2001). 
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which is consistent with the upper limit on T from the b— value obtained by Sembach et al. (2001). 
Figure 4 also shows that the absence of detectable Si IV and C IV absorption in this system is fully 
consistent with the model. At log U ~ —3.35, the predicted Si IV and C IV columns are more 
than an order of magnitude below the observational limits. The model is consistent with the other 
upper limits as well, including the upper limit on N(C III) set by Sembach et al. (2001). We note 
that the model is required to produce the observed H I column density, log N(R I) = 15.85, and 
consequently as U is increased and the gas becomes more highly ionized, the thickness of the slab 
is increased to match the observed H I column. This is why iV(Si II) and N(Si III) both increase 
with increasing U even though the ion fraction of Si II decreases with increasing U over the range 
shown in Figure 4. 

It is interesting to note that silicon, an a— group element, is overabundant compared to C 
by ~ 0.4 dex in this model. We initially attempted to reproduce the observed column densities 
assuming that the relative heavy element abundances in the Virgo absorber follow the solar pattern 
as well (i. e., [Si/C] = 0.0), but we found that this could not be done. This is apparent from a brief 
inspection of Figure 4: when the Si Ill/Si II ratio is reproduced, the model produces too much C II 
if [Si/C] = 0.0 (multiphase models presented below require a somewhat lower Si overabundance, 
but the best fit in these models still yields [Si/C] ~ 0.2). We discuss the implications of this result 
in §4. 

As shown by Weymann et al. (1995) with regard to H ionization, local sources of photoion- 
ization, such as X-rays from the hot gas in the Virgo cluster or UV flux escaping from nearby 
galaxies, are unlikely to be important compared to photoionization by the background UV flux 
from QSOs. We confirm that the heavy element results presented above are insensitive to these 
issues. For example, if we include X-ray flux from the hot gas (following §3.1.4 in Weymann et al. 
1995) in the CLOUDY model, we find only miniscule changes in the column densities of interest 
compared to Figure 4. Similarly, if we assume that the radiation field is dominated by the UV flux 
from stars, we obtain very similar constraints on the heavy element abundances. Even lower Si IV 
and C IV column densities are predicted in this case, so the observational constraints in Table 1 are 
easily satisfied. However, if stellar flux were competitive with the QSO background, the gas density 
would be higher; the combined flux from the QSO/AGN background and the stellar sources leads 
to a higher photon density, so the gas density must also be higher to produce the same ionization 
parameter. 

We next consider the possibility that the absorption occurs in a multiphase medium in which 
the C II, Si II, and Si III lines do not entirely come from the same gas. As shown in Figure 5, the 
centroids of these lines (see also column 5 in Table 1) marginally support this hypothesis: the C II 
and Si II lines are blueshifted by ~5 km s" 1 compared to the centroid of the Si III absorption. Also, 
the Si III profile appears to show two components with the weaker component more closely aligned 
with the Si II and C II lines. However, the Si II and C II lines are very weak and therefore could 
be misleading due to noise. Consequently, we conclude that while there are reasonable indications 
that this is a multiphase absorption system, this requires confirmation with higher S/N data. 
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Fig. 4. — Heavy element column densities predicted for a Lya absorber, photoionized by the UV 
background from QSOs (from Haardt & Madau 1996) with log iV(H I) = 15.85 and heavy element 
abundances indicated at upper left, compared to the observed column densities at z a b s = 0.00530. 
The model column densities are indicated with solid lines with small symbols while the observed 
columns are shown with larger symbols with la error bars (or arrows in the case of 3a upper limits). 
The species corresponding to the symbols are shown in the key at the right. 
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Even if the absorber is a multiphase medium, we can still derive useful estimates of the heavy 
element abundances because the H I, C II, and Si II absorption should predominantly originate in 
the same, lower-ionization, phase. To show this, we plot in Figure 6 the ion fractions of H I, C II, 
Si II, and Si III vs. log U from the same model shown in Figure 4. We also show the gas temperature 
predicted by the model over this range of U. From this figure one can see that the ion fractions of 
H I, C II, and Si II show very similar trends as the gas becomes more ionized, and it is reasonable 
to assume that the lines of these three species arise in the same gas. Furthermore, Si III does not 
follow the lower ionization species, and it would not be surprising to find that there is somewhat 
more ionized gas that produces a Si III line that is substantially stronger than the Si II line. The 
bulk of the H I detected in higher Lyman series lines should be affiliated with the Si II/C II gas. 
However, the higher-ionization Si Ill-bearing gas may produce weaker H I absorption that may be 
detectable in the Lya transition (unless it is confined to the strongly saturated portion of the Lya 
profile). More importantly, it is evident from Figure 6 that some Si III arising in the Si 11/ C II 
phase is likely to be detectable, and this may explain the indications of a weaker component in the 
Si III profile at the velocity of the Si II and C II lines. The column density of this weaker component 
is poorly constrained due to blending of the two components and the S/N of the data, but fitting 
two components to the Si III profile with the Fitzpatrick & Spitzer (1997) profile-fitting code and 
the line spread functions from the STIS Handbook (Leitherer et al. 2001) yields log iV(Si III) 
pa 11.9 ± 0.2 for the weaker component. 

Assuming the H I, C II, and Si II lines go together and adopting log N(Si III) « 11.9 ± 0.2 for 
the Si 11/ C II-bearing gas, we find from the photoionization models that the following abundances 
are consistent with the observed column densities within the la uncertainties: 

[c/h] = -1.218:1, (i) 

and 

[Si/C] = +0.2 ±0.1, (2) 

with 

log n H = -2.8 ± 0.3. (3) 

At the current S/N level, the presence of a more highly ionized phase that produces an additional 
Si III component (discussed in the previous paragraph) does not conflict with the upper limits on 
Si IV and C IV. However, at higher S/N levels, C IV and/or C II absorption associated with this 
more highly ionized component should be detectable. It is possible that some of the C II absorption 
that we have already detected is associated with the more highly ionized extra Si III component. 
If so, then we have overestimated [C/H] and we have underestimated [Si/C]. Decreasing [C/H] and 
increasing [Si/C] to compensate for this would only strengthen our interpretation of the nature of 
this absorber (§4). 

The model has some interesting implications. With these constraints on ran, the model predicts 
11,000 <T< 14,000 K (see Figure 6), and the cloud pressure P/k « 40 K cuT 3 . However, the 
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Fig. 5. — Continuum-normalized absorption profiles of the C II, Si II, and Si III lines at z a b s = 
0.00530, plotted on an expanded velocity scale to show the apparent differences in the line centroids 
and component structure. For reference, the tick mark indicates a velocity of —9 km s" 1 . 
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Fig. 6. — Logarithmic ion fraction (e.g. Si Il/Sitotai) from the model shown in Figure 4 vs. log U 
(bottom axis) and log «h (top axis) for H I (stars), C II (filled squares), Si II (filled triangles), and 
Si III (open triangles). The temperature of the gas is also shown with a dashed line using the linear 
scale on the right-hand axis. 
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thickness of the gas slab is only ~70 pc for log nu = —2.8. This is vastly smaller than expected 
for a gravitationally confined Lya cloud with an ordinary ratio of baryons to dark matter; a 
photoionized, self-gravitating cloud with log iV(H I) = 15.85, the inferred temperature, and a gas 
fraction f g pa Qb/Q m would be several orders of magnitude larger (see Schaye 2001). The inferred 
absorber thickness is somewhat sensitive to the ionization parameter. For example, if we decrease 
the density to log ran = —3.1 (i.e., log U = —3.35), then the thickness increases to ~ 350 pc. This 
is still much smaller than expected for a self-gravitating cloud. The cloud mass is uncertain due to 
the unknown geometry and impact parameter. If spherical, the cloud has a rather small baryonic 
mass on the order of 10 — 400 M Q , but the mass could be substantially larger with other geometries. 
Similar properties have recently been derived for a subset of weak Mg II systems at z ~ 1 (Rigby, 
Charlton, & Churchill 2002). 

The small line-of-sight size implied by the model indicates that either (1) the fraction of the 
absorber mass in the form of baryonic gas is extremely small (f g ~ 10~ 6 — 10~ 5 ), (2) the absorber 
is pressure-confined by an external medium, or (3) the gas is out of equilibrium and will rapidly 
expand or evaporate. A small / g could arise in a small dark matter halo that expelled most of 
its gas in an early episode of star formation. Mac Low & Ferrara (1999) have shown that SN can 
expel most of the ISM from a very low-mass object. Alternatively, the baryons may have been 
photoevaporated from a small dark matter halo during the reionization epoch (e.g., Klypin et al. 
1999; Barkana & Loeb 1999). Assuming the current gas mass is ~ 10 M Q , the required f g would 
imply that the mass of the dark matter halo is 10 6 — 10 7 M . A substantial fraction of the gas 
can be removed from a dark matter halo in this mass range (see references above). However, the 
baryon removal process must be remarkably efficient, and the gas which remains in the dark matter 
halo must be in a relatively quiescent state (there is some evidence of component structure in the 
absorption lines, but the profiles are quite simple compared to many QSO absorption systems). 
Furthermore, the fact that the baryonic mass (and therefore the implied dark matter halo mass) 
could be substantially higher with other (perhaps more likely) geometries such as a filament or 
sheet makes baryon evacuation more difficult and the required f g value more improbable. 

It seems more likely that the absorber is pressure confined by an external medium. For purposes 
of illustration, Fi gure 7 shows how the Lyo: profile at Zabs 

= 0.00530 can easily accomodate a broad 
component with a 6— value as large as 41 km s — 1 . If due to the external confining medium, this 
broad-component Doppler parameter implies that T ext < 1 x 10 5 K, and therefore the external 
medium would need a density n cxt > 4 x 10~ 4 cm -3 to confine the metal-bearing gas with the 
pressure derived above. This is a plausible density for a galaxy halo or the intracluster medium 
in a galaxy group. However, the fit shown in Figure 7 is not unique, and it is possible that a 
substantially broader H I component is present. A sufficiently broad and weak component could 
be hidden in the noise, particularly since the strength of such a feature could be inadvertently 
reduced by the continuum fitting process. Consequently, T cxt could be substantially higher and 
n cxt accordingly lower. 
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Fig. 7. — Broad-component fit to the Lya line at z a bs = 0.00530. In both panels, the observed 
Lya profile is shown with a histogram plotted versus rest-frame velocity (v — km s ^at z&\ys 
= 0.00530). In panel (a), a broad component with b = 41 km s _1 and log iV(H I) = 13.73 is 
overplotted on the observed profile with a thick solid line, and the main component with b = 16.1 
km s _1 and log iV(H I) = 15.85 (as derived by Sembach et al. 2001 from a curve-of-growth analysis) 
is shown with a thick dashed line. Panel (b) shows a two-component fit to the Lya line including 
both components shown in panel (a). 



-19- 



3.2. z ahs = 0.00337 

We now constrain the metallicity of the 3C 273 Lya cloud at z a b s = 0.00337, first assuming 
the gas is collisionally ionized (§3.2.1), and then considering the possibility that it is photoionized 
(§3.2.2). The only detected species are H I and O VI. From a curve of growth analysis, Sembach 
et al. (2001) estimate that this system has log iV(H I) = 14.41±0.10 at this redshift. However, the 
absence of C IV provides a useful constraint for ionization models, and even with only two detected 
species, we can still place limits on the metallicity of the gas as follows. The logarithmic oxygen 
abundance is given by 

"0 
H 

where / is the ion fraction and (O/H) is the solar oxygen abundance (we take log (0/H)q = —3.26 
from Holweger 2001). In general, it is usually possible to minimize /(H l)//(0 VI), and this sets 
a lower limit on [O/H] since the other quantities in eqn. 4 are measured. In practice, it can be 
difficult to ascertain the 7V(H I) that should be assigned to the O Vl-bearing gas; there is evidence 
that some O VI absorbers are multiphase entities (e.g., Tripp, Savage, & Jenkins 2000), and some 
of the H I absorption may arise in a lower-ionization phase which does not contain O VI. 



3.2.1. Collisional Ionization 

We first consider collisional ionization. If the O Vl-bearing gas is collisionally ionized, then the 
lower limit on O Vl/C IV at z a b s = 0.00337 provides a useful lower limit on the gas temperature. 
Figure 8 shows the C IV and O VI ion fractions in collisional ionization equilibrium (Sutherland &: 
Dopita 1993) and in nonequilibrium radiatively cooling gas (Shapiro & Moore 1976) for T > 10 5 K. 
We see that at T < 10 5 2 K, the ionization of the gas strongly favors C IV if the gas is in equilibrium; 
in this case the C IV ion fraction is one to several orders of magnitude higher than the O VI ion 
fraction. To set a lower limit on T using the lower limit on O VI/C IV observed at z a b s = 0.00337, 
we must make a further assumption about the intrinsic (O/C) abundance in this absorber. In 
general, we use the solar relative abundances of Holweger (2001). However, we presented evidence 
in the previous section that a— elements are overabundant with respect to carbon in at least one 
Virgo Lya cloud. Since O usually follows the a— element abundance patterns (McWilliam 1997), 
we will assume [O/C] = 0.3; this will set a more conservative lower limit on T. In this way, we 
derive T > 10 5 29 K from the O VI/C IV limit for the equilibrium case. C IV is less strongly favored 
at lower temperatures in the nonequilibrium model, but nevertheless the O VI/C IV ratio requires 
a similar lower limit on T; from the nonequilibrium model we obtain T > 10 518 . 17 



17 We note that in the nonequilibrium radiatively cooling gas calculations of Schmutzler & Tscharnuter (1993), the 
C IV ion fractions are substantially lower, and the O VI/C IV constraint can be satisfied in much cooler gas. The 
difference is evidently due to atomic data since a comparison of the equilibrium calculations of Shapiro & Moore 
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Fig. 8. — Ion fractions of C IV (squares) and O VI (stars) as a function of temperature for gas 
in collisional ionization equilibrium (solid lines, from Sutherland k, Dopita 1993) and nonequilib- 
rium raditiatively cooling gas (dashed lines, from Shapiro k, Moore 1976). In the nonequilibrium 
calculation, the gas is initially in equilibrium at 10 6 K and then cools isochorically. 
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To constrain the H I column associated with the O VI gas, we turn to the Lya profile. If the 
O VI is due to hot, collisionally ionized gas, then the corresponding H I absorption line should 
be broad. We note that since the observed O VI/C IV limit requires a similar temperature in the 
equilibrium and non-equilibrium radiatively cooling gas models, in this absorber a broad component 
is expected even if the gas was suddenly shock- heated and subsequently cooled faster than it could 
recombine (as modeled by Shapiro & Moore 1976). The H I b— value derived by Sembach et al. 
(2001) from the COG analysis, b = 30lg'g km s^ 1 , is not compatible with the limits on T from 
O VI/C IV; this b implies a temperature upper limit (T < 55,000 K) at which C IV should be 
substantially stronger than O VI if collisionally ionized. Even at the 3er upper limit, this 6— value 
implies that C IV should be stronger than or at least comparable to the O VI (see Figure 8). 

However, inspection of the top panel of Figure 3 reveals evidence of structure in the Lya profile. 
It is apparent that a sufficiently broad H I component associated with the O VI is allowed by the 
data as long as the Lya profile probes multiple phases of the gas. To set a lower limit on [O/H] at 
Zabs = 0.00337, we must minimize /(H l)//(0 VI). In collisional ionization equilibrium, this occurs 
at log T » 5.5, where log /(H I)//(0 VI) = -5.14 (Sutherland & Dopita 1993). Figure 9 shows 
a fit to the Lya profile including a broad component consistent with log T = 5.5 (which requires 
b « 72 km s" 1 ). The upper panel in Figure 9 shows the broad component only, with b = 72 km 
s _1 and log iV(H I) = 13.42. To obtain a good fit to the overall profile, an additional, narrower H I 
component due to a separate, cooler phase must be added, and the lower panel in Figure 9 shows 
a fit including the broad component in Figure 9a as well as a cooler component with 6 coo i = 33 
km s _1 and log iV coo i(H I) = 14.20. The two-component model shown provides an acceptable fit to 
the Lya profile and is compatible with the Ly/3 profile recorded with FUSE. So, with log iV(H I) 
= 13.42, log /(H l)//(0 VI) > -5.14, and log N(0 VI) = 13.32 for the hot, O VI gas, we obtain 
[O/H] > -2.0. 

Figure 8 compares the equilibrium O VI and C IV ion fractions to those from a particular 
non-equilibrium model. Several other non-equilibrium collisionally ionized gas models have been 
developed such as conductive interfaces or turbulent mixing layers between very hot gas and cool 
clouds (see §8 in Sembach et al. 1997 for a brief review). Some of these other models, e.g., the 
conductive interface model of Borkowski et al. (1990), can satisfy the O VI/C IV limit in the z a b s 
= 0.00337 Virgo absorber. However, we note that the turbulent mixing layer model of Slavin et 
al. (1993) does not produce enough O VI; all of the parameter combinations considered by Slavin 
et al. lead to O VI/C IV < 1. Moreover, the O VI column densities from a single turbulent mixing 
layer are more than an order of magnitude smaller than the observed O VI column in this system, 
so a large number of layers would be needed to explain the observations. 



(1976) to Schmutzler & Tscharnuter shows the same C IV discrepancy. We favor Shapiro & Moore since the C IV 
ion fractions from the large majority of collisional ionization equilibrium calculations (e.g., Arnaud & Rothenflug 
1985; Shull & van Steenberg 1982; Sutherland & Dopita 1993; Mazzotta et al. 1998) agree with the Shapiro & Moore 
steady state results. 
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Fig. 9. — Broad-component fit to the Lya line at z a \> s = 0.00337. As in Figure 7, the observed 
Lya profile is shown with a histogram plotted versus rest-frame velocity (v — km s ^at 
= 0.00337). In panel (a), a broad component with b = 72 km s _1 and log iV(H I) = 13.42 is 
overplotted on the observed profile. Panel (b) shows a two-component fit to the Lya line including 
the broad component shown in panel (a). 
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3.2.2. Photoionization 

We next derive constraints on the metallicity of the z a b s = 0.00337 system assuming the O VI 
is photoionized. In this case, we assume that all of the neutral hydrogen derived from the curve-of- 
growth analysis (log 7V(H I) = 14.41) is associated with the O VI phase since this provides the most 
conservative lower limit on [O/H]. In this circumstance, the O Vl-bearing gas is consistent with the 
narrow core of the Lya line because it can be much cooler than in the collisionally ionized case. We 
use the photoionization models from §3.1 to determine the minimum value of /(H l)//(0 VI), but 
we impose the additional stability requirement that the size of the absorbing cloud not exceed the 
Jeans length. 18 With these assumptions, we find log /(H l)//(0 VI) > —4.12 with log nn > —5.5 
and absorber thickness L < 1.5 Mpc. Therefore we obtain a similar lower limit to the oxygen 
abundance, [O/H] > —2.0. A smaller L requires a higher metallicity; L < 450 kpc requires [O/H] 
> —1.5, for example. In principle, the absence of C IV provides an upper bound on the metallicity, 
but this is compromised by the uncertain iV(H I) associated with the O VI: if the absorber is 
multiphase and log iV(H l)ovi < 14.41, then this increases the lower limit on [O/H] (see eqn. 4). 
Consequently, we can place a firm lower limit on the oxygen abundance but we cannot place a 
strong upper limit on [O/H]. 

4. Discussion 

We now offer some comments on the implications of these abundance measurements. We first 
discuss the z^bs = 0.00530 absorber. Using the abundance patterns observed in the Milky Way 
ISM as a reference, we consider whether the Si overabundance in this system could be due to 
depletion by dust (§4.1) or an intrinsic a— element overabundance (§4.2); we favor the latter. We 
then discuss the environment in which the absorber is found and the nature of the absorption 
system (§4.3). Finally, we make some brief comments on the detection of O VI at the redshift of 
the Virgo supercluster (§4.4). 

4.1. Depletion of C and Si by Dust 

We have obtained gas-phase abundance measurements for carbon and silicon at z a b s = 0.00530. 
In the ISM of the Milky Way, both of these elements can be depleted by incorporation into dust 
grains. Therefore the total abundances, i.e., gas-phase + dust-phase, could be higher than the 
numbers reported in §3. In the Galactic ISM, Si is usually more substantially depleted than C, 
which would leave little room for dust depletion in the Virgo absorber studied here. For example, 
toward C Oph, [C/H] = -0.4 while [Si/H] = -1.3 (e.g., Sofia, Cardelli, k Savage 1994; Snow & 
Witt 1996). Before proceeding with this comparison, we note an important caveat. In §3, we 



This stability requirement is easily satisfied in the best-fitting models in §3.1. 
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modeled the ionization of the Virgo Lya absorber in order to derive the Si and C abundances in 
the gas. In the Milky Way ISM, such modeling is usually not applied. Instead, measurements of 
the column densities of ionization stages that are dominant in H I regions are taken to provide a 
good estimate of the elemental abundance, e.g., [C/H] = log N(C II)/JV(H I) — log (C/H) Q . In this 
way, gas-phase abundances of species such as Si, Cr, and Fe are found to be substantially reduced 
compared to abundances of elements such as O, S, and Zn, species that are depleted lightly or not 
at all. This is usually assumed to be primarily due to depletion by dust. However, in some cases 
ionization effects may also play a role (e.g., Jenkins et al. 1998; Sofia & Jenkins 1998; Jenkins et 
al. 2000). In some cases, neutrals are included in the abundance estimate (e.g., both C I and C II 
are measured for the purpose of estimating [C/H]), but a larger error can result from the fact that 
iV(H I) can underestimate the total hydrogen density associated with a particular heavy element 
because the dominant ionization stage can exist in gas where the H is largely ionized. Sembach et 
al. (2000) have modeled the ionization of the "warm ionized medium", but ionization corrections 
are often neglected in the ISM literature. 19 With this in mind, we note that in interstellar gas in 
the disk, [Si/H] ranges from —0.35 to —1.31 (Sembach & Savage 1996) while Sofia et al. (1998) 
find [C/H] ss —0.45 over a wide range of physical conditions. In the Galactic halo, Si depletion is 
reduced, e.g., Sembach & Savage (1996) find that [Si/H] ranges from —0.09 to —0.47 on halo sight 
lines. Depletion of carbon in halo gas has not been extensively studied. 

The Galactic carbon and silicon abundances reported by Sofia et al. (1998) and Sembach 
& Savage (1996), respectively, are not derived from the same interstellar clouds, and it is not 
entirely clear how to combine the two samples to examine the [Si/C] abundance patterns that 
might result from depletion by dust. There are concerns about estimating [Si/C] trends from Si 
and C measurements in different clouds. For example, if dust destruction liberates silicon leading 
to the generally larger Si abundances in halo gas (see above) , the same processes could increase the 
C abundances in the halo as well. Also, the overall metallicity of the ISM might not be constant, 
in which case the [Si/C] ranges derived by combining the results of Sofia et al. (1998) and Sembach 
& Savage (1996) could be misleading. To circumvent these issues, we have assembled in Table 3 a 
sample of interstellar clouds in which both Si II and C II have been reliably measured in the same 
cloud, including new measurements for two sight lines (see Appendix). We have limited this sample 
to measurements made with the highest resolution modes of GHRS and STIS (i.e., FWHM < 3.5 
km s _1 ) in order to minimize potential confusion due to unresolved saturation. From this table, 
we see that the majority of the clouds in this sample have [Si/C] <C 0.0. In four cases in Table 3, 
[Si/C] is consistent with the solar ratio or somewhat larger. However, in two of these four cases, 
the C II line approaches zero intensity in the line core, and N(C II) could be underestimated due to 
unresolved saturation, as noted by Lehner et al. (2001) [all of the Si II measurements are unlikely 
to be underestimated due to saturation]. Increasing N(C II) would decrease [Si/C]. The remaining 
two cases with high [Si/C] values (both from the DI 1388 sight line) are evidently uncertain; for 



Ionization corrections are most substantial in low density clouds; in many Milky Way clouds the densities are 
high enough so that the corrections are likely to be small. 



-25- 



these clouds Lehner et al. report rough measurements without error bars. It is also worth noting 
that there are tentative indications of a— element overabundances in some high-velocity clouds 
(e.g., Richter et al. 2001), and the high values of [Si/C] in the HVCs toward DI 1388 may be due 
to intrinsic Si overabundances. We conclude that the sample in Table 3 presents no compelling 
evidence that dust depletion could lead to a gas-phase silicon overabundance with respect to carbon. 
On the contrary, it seems likely that dust would produce an underabundance of Si in the gas phase. 

4.2. a— Element Abundances 

It is quite possible that silicon is intrinsically overabundant in the Virgo Lya cloud at z a b s 
= 0.00530, and that dust depletion affects the abundances by only a small amount if at all. This 
interpretation has interesting implications. It is well known that the a— group elements, including 
silicon, are overabundant relative to iron by roughly 0.3 dex in low metallicity stars in various 
locations in the Milky Way (Wheeler et al. 1989; McWilliam 1997, and references therein). Fol- 
lowing Tinsley (1979), this is generally believed to indicate that the heavy element enrichment is 
dominated by the output from Type II supernovae (SN II). The nucleosynthetic origins of carbon 
are somewhat complex, and carbon abundance trends in stars are highly uncertain, but McWilliam 
(1997) notes that there are no compelling indications that C is overabundant (relative to Fe) in 
low-metallicity Galactic halo stars. Furthermore, some galactic chemical evolution models can pro- 
duce overabundances of Si relative to C, with [Si/C] as high as ~0.5 (e.g., Timmes et al. 1995). 
Therefore there is reason to expect [Si/C] ~ 0.3 in low-metallicity objects such as this Virgo Lya 
cloud. 



4.3. Nature of the Absorber at z a b s — 0.00530 

There is also long-standing evidence that a— group elements are overabundant (with respect 
to iron) in the X-ray emitting hot intracluster medium (ICM) in several clusters including Virgo 
(e.g., Canizares et al. 1982; Mushotzky et al. 1996; Matsumoto et al. 1996). This has usually 
been ascribed to galactic winds driven predominantly by SN II (e.g., Loewenstein &: Mushotzky 
1996, but see the discussion of caveats in Gibson, Loewenstein, & Mushotzky 1997). More recent 
analyses have concluded that the Fe abundances increase in the central regions of several clusters 
while the a— group abundances remain more-or-less constant (e.g., Finoguenov, David, k, Ponman 
2000; Kaastra et al. 2001). This likely indicates that ejecta from Type la SN contribute to the ICM 
enrichment in the central regions where the cD and other old elliptical galaxies are located, but SN 
II ejecta from early "protogalactic" winds dominate the enrichment in the rest of the ICM. The 
sight line to 3C 273 is > 10.3° away (i.e., projected distance > 3 Mpc) from M87 and the X-ray 
emitting ICM in the Virgo cluster. Consequently, it is interesting that the Virgo Lya absorber 
at z a bs = 0.00530 also shows evidence of an overabundance of one of the a elements. This may 
indicate that the early SN II winds affect a very broad region of the cluster, a region much more 
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Table 3. High-Resolution [Si/C] measurements in the Milky Way ISM 
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b All [Si/C] values from the literature have been recalculated assuming the Holweger (2001) solar abun- 
dances. Error bars include uncertainties in the solar reference abundances reported by Holweger (2001) as 
well as the uncertainties in the column density measurements. 

c Comment: L = low density sight line within the "Local Bubble", D = high density sight line, F = 
forbidden velocity cloud (i.e., velocity is outside of the expected range for normal Galactic rotation for this 
direction and background star distance), O = star is behind "Orion's Cloak" (Cowie, Songaila, & York 
1979), U = due to the substantial distance of the background star, the location of the absorption system 
is uncertain. Based on the observed velocities of Magellanic Bridge gas, these absorbers are more likely to 
be associated with the Milky Way than the Magellanic Bridge (Lehner et al. 2001). 

d Column density measurements from: (1) Cry & Jenkins (2001), (2) Cardelli et al. (1993), (3) Cardelli 
ct al. (1994), (4) Welty et al. (1999), (5) Howk, Savage, & Fabian (1999), (6) this paper (see Appendix), 
(7) Lehner, Keenan, & Sembach (2001). 

°The C II absorption approaches the zero flux level in the core of the line, and N(C II) could be 
underestimated due to saturation. Consequently, these measurements are most conservatively treated as 
upper limits. The upper limits are the nominal best values (and corresponding la uncertainties) from the 
original sources. 

'This star is located within the gaseous "Magellanic Bridge" between the Large Magellanic Cloud and 
the Small Magellanic Cloud. 

8 Column density uncertainties are not reported in the original reference. 
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extended than the easily observed hot gas. This also has implications for the issues briefly noted 
in §1. For example, if the intergalactic gas within Virgo is enriched by dynamical processes such as 
ram-pressure stripping, the abundance patterns should reflect a significant contribution from SN la 
(Gunn & Gott 1972). The observation of SN II abundance patterns favors supernova driven winds 
as the source of enrichment. Furthermore, if the z a b s = 0.00530 absorber is a small blob entrained 
in the hot gas of a galactic wind, then the hotter wind gas could provide the pressure required to 
confine the absorber, as discussed in §3.1. 

We note that analogous abundance trends have been observed in other contexts in low- and 
high-redshift absorption systems. For example, Richter et al. (2001) have measured abundances in 
the high- velocity cloud (HVC) Complex C near the Milky Way using absorption lines detected in 
the spectrum of PG1259+593 and a well-constrained COG. They report [O/H] = -l.Olojj, [Si/H] 
= —0.9 ±0.3, and [Fe/H] = -1.3^. It is interesting to note the overabundance (albeit marginal) 
of the a— group elements O and Si with respect to Fe in this object which, based on its large 
angular extent, is probably relatively close to the Galaxy (d ~ 10 kpc; Wakker et al. 1999). In this 
case the abundance trends could be entirely due to dust, which depletes iron more strongly than 
O and Si and can thereby mimic a overabundances. The Virgo absorber at z a ^ s = 0.00530 could 
be analogous to this HVC, although its H I column density is substantially lower. 

There is also evidence of a— group overabundances in high-redshift damped Lya absorption 
systems (e.g., Lu et al. 1995, 1996; Prochaska & Wolfe 1999). In comparisons of this Virgo Lya 
absorber to higher redshift systems, it is important to note that low— z and high— z absorbers with 
the same iV(H I) are probably not dynamically equivalent. In cosmological simulations, a Lya 
system with a given column density at z ~ has the same overdensity as a system at z ~ 3 with 
a larger iV(H I) by a factor of ~ 50 (Dave et al. 1999). Therefore this Virgo Lya absorber may be 
more analogous to a Lyman limit system at high redshift. 

However, damped Lya and Lyman limit absorbers are generally thought to be associated with 
individual galaxies. Despite its location within the Virgo supercluster (or perhaps because of this), 
it is difficult to attribute the z a b s = 0.00530 system to such a galaxy. After the Virgo Lya absorbers 
in the 3C 273 spectrum were discovered (Morris et al. 1991; Bahcall et al. 1991), several studies 
of the relationships between these systems and galaxies were carried out making use of galaxy 
redshift measurements (Salzer 1992; Morris et al. 1993; Hoffman, Lewis, & Salpeter 1995; Salpeter 
& Hoffman 1995; Bowen, Blades, & Pettini 1996; Grogin & Geller 1998; Impey, Petry, k Flint 1999; 
Penton, Stocke, & Shull 2000, 2002), deep optical imaging (Morris et al. 1993; Rauch, Weymann, 
& Morris 1996) and 21cm imaging (van Gorkom et al. 1993). Despite these intensive searches 
with good sensitivity to faint and low surface-brightness galaxies, no galaxies have been found with 
impact parameter p less than 100 kpc near the redshift of the z a b s = 0.00530 absorber. Damped 
Lya and Lyman limit absorbers are generally believed to have cross sections significantly smaller 
than 100 kpc (e.g., Steidel 1993). 

There are some known galaxies near z = 0.00530 with p 3> 100 kpc (see Table 1 of Salpeter & 
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Hoffman 1995 as well as Grogin, Geller, & Huchra 1998). For example, as noted by Salzer (1992), 
UGC 7549 (NGC 4420) is within 100 km s" 1 of this absorption system with p = 270 kpc. UGC 
7549 is notable only because it is the brightest galaxy in the vicinity of 3C 273; in other respects it 
appears to be an ordinary spiral galaxy However, there are at least two other galaxies at similar 
projected distances and velocity separations: UGC 7612 (p = 244 kpc, | At; |= 80 km s _1 ) and 
the low surface brighness galaxy UGC 7642 (p = 241 kpc, | Av \ = 50 km s" 1 ). H I 1225+01 is 
at a somewhat larger, but still plausible, velocity separation (p = 167 kpc, | Av |= 300 km s _1 ). 
H I 1225+01 is a large cloud of H I detected in 21cm emission (Giovanelli & Haynes 1989), part 
of which envelops a small dwarf galaxy (Salzer et al. 1991). This object is intriguing because the 
H I envelope has a complex, two-component morphology, and the dwarf galaxy shows regions of 
active star formation (Giovanelli, Williams, & Haynes 1991). This is the type of galaxy that may 
be able drive a wind with which the z a b s = 0.00530 system could be associated (see Heckman et al. 
2001 for an observational example of such a wind). 20 Theuns, Mo, & Schaye (2001) have suggested 
that dwarf galaxy winds "may break up into clouds, which coast to large distances [~300 kpc] 
without sweeping up a significant fraction of the IGM." The small size and mass derived for the 
-Zabs = 0.00530 Lya cloud is qualitatively consistent with this suggestion, and as noted above, a 
small cloud entrained in a wind could be naturally pressure confined. However, while it is tempting 
to assign the gas in this absorber to one of these galaxies, we see no clear method to determine 
which one, if any, is the correct assignment. Given the complex velocity field within the Virgo 
cluster, it remains possible that some of these galaxies are interlopers, i.e., the three-dimensional 
distances from the Lya cloud to the galaxies are larger than implied by their velocity differences. 
Furthermore, there is evidence that the two 3C 273 absorbers in Virgo and their nearest galaxies 
are part of a large-scale filament (Penton et al. 2001), in which case the absorption line system 
may not be associated with any particular galaxy but rather is part of the general intragroup gas. 

Lacking a specific galaxy to which this absorption system can clearly be assigned, it is interest- 
ing to consider whether this Virgo absorber might be a relatively pristine gas cloud that has only 
been enriched by a putative first generation of stars, i.e., the "Population III" stars (e.g., Qian, 
Sargent, &: Wasserburg 2002). Using the Population III yields calculated by Heger & Woosley 
(2002), Oh et al. (2001) have argued that explosions of the first very massive stars could produce 
overabundances of a— group elements with respect to carbon, exactly as observed in the Virgo Lya 
cloud. However, Oh et al. (2001) have also shown that while there is evidence of Pop III enrich- 
ment in low-metallicity stars with Z < 10~ 3 Z Q , at higher metallicities enrichment from normal 
stars appears to be dominant. Therefore the metallicity derived for the Virgo absorber (see §3) 
appears to be too high to be consistent with this hypothesis. 



20 However, we note that mechanisms have been proposed (e.g., Binney 2001) which would drive winds from luminous 
larger galaxies such as UGC 7549. 
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4.4. Warm/Hot Gas in Virgo 

In the process of constraining the oxygen abundance in the Virgo O Vl/Lya absorber at z a bs = 
0.00337, we have explored the physical conditions of this system (§3.2). We now briefly comment on 
these physical condition constraints. Based on hydrodynamic simulations of cosmological structure 
growth, it has been suggested that at the present epoch, 30-50% of the baryons are in low-density, 
shocked gas at 10 5 — 10 7 K, the so-called "warm/hot intergalactic medium" (Cen & Ostriker 1999a; 
Dave et al. 2001). Lithium-like O VI is useful for probing such gas since its ion fraction is maximized 
at T « 10 5 ' 5 K in collisionally ionized gas (see Figure 8), and low— z O VI absorbers do harbor a 
substantial quantity of baryons (Tripp, Savage, & Jenkins 2000; Savage et al. 2002). However, the 
nature of O VI absorbers is not yet clear. Are these systems collisionally ionized or photoionized? 
Do they arise in unvirialized, large-scale filaments (as predicted by the cosmological simulations, 
see Cen et al. 2001; Fang k, Bryan 2001), higher density galaxy groups (Mulchaey et al. 1996), 
individual galaxies, or galactic winds? 

The 3C 273 O VI-Lya absorber at 

^abs — 0.00337 provides some interesting clues. The absorp- 
tion line properties are consistent with expectations for collisionally ionized, warm-hot intergalactic 
gas at T > 10 5 K. As Lya systems go, this absorber is located in a region of relatively high galaxy 
density (Grogin & Geller 1998). Given the low O VI equivalent width, this absorber may not be 
located in the typical low- overdensity environment expected based on the cosmological simula- 
tions (see, e.g., Figure 4 in Cen et al. 2001), but this is only one case, and it is quite likely that 
analogous O VI systems can be found in the simulations. Like the z a b s = 0.00530 system, there 
are no galaxies very near the sight line at z a b s = 0.00337 (Hoffman et al. 1998, and references 
therein), and it is not clear whether or not this absorber is in a virialized region. Additional galaxy 
redshift measurements currently underway (K. McLin et al., in preparation) will help illuminate 
the nature of this absorption system. Other O VI absorbers at higher redshifts have line properties 
similar to the z a b s = 0.00337 system, e.g., similar lower limits on O VI/C IV and evidence of broad 
components in the Lya profiles (e.g., Tripp & Savage 2000; Tripp et al. 2001; Savage et al. 2002). 
It is interesting to note that in all of these absorbers, the broad H I lines are juxtaposed with at 
least one narrow H I line. This might be expected if the broad-component gas was shock-heated 
as high-velocity gas collided with the lower-ionization gas that gives rise to the observed narrow 
component. On the other hand, if the broad H I component is due to Hubble broadening over a 
long pathlength (as may be required if the O VI is photoionized, see Tripp et al. 2001), then there 
is no particular reason to expect the observed correspondence between the broad and narrow H I 
lines. The 3C 273 O VI system at -^abs — 0.00337 is particularly valuable for understanding these 
absorbers because the low redshift allows an intensive investigation of the absorber environment, 
and because the QSO is the brightest known. 
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5. Summary 

We have obtained good constraints on the abundances of carbon and silicon in a Lya absorption 
system at z a b s = 0.00530 in the spectrum of 3C 273; this absorber is located within the southern 
extension of the Virgo cluster. Based on ionization modeling, we find that the carbon abundance 
in the gas phase is roughly 1/15 of the solar abundance, and Si, an a— element, is overabundant 
relative to C by ~0.2 dex. The H density is ~ 0.0015 cm -3 , the absorber thickness is ~70 pc, and 
the gas pressure p/k ~ 40 cm~ 3 K. The small line-of-sight thickness suggests that the absorber is 
pressure confined by an external medium or is out of equilbrium and will soon expand or evaporate 
(a self-gravitating cloud would require a remarkably high dark matter to baryon ratio). The Si 
overabundance is probably not due to dust depletion. Instead, this likely reflects an intrinsic 
abundance pattern. Overabundances of a— group elements have been observed in the hot X-ray 
emitting gas in the Virgo cluster, and these abundance patterns are believed to be produced by 
galactic winds containing predominantly SN II ejecta. We favor such an explanation for this 3C 
273 absorber as well because there are no galaxies particularly nearby, as might be expected in 
other scenarios, but there are galaxies at projected distances of ~200 kpc which could plausibly 
drive an escaping wind. The projected distance from the 3C 273 sight line to the region detected in 
X-rays is ~3 Mpc, and therefore the putative SN II winds may be widespread in the Virgo cluster. 
If this result is found for other sight lines, SN II winds, perhaps "protogalactic" winds, may not 
be unique to the central regions of clusters and may in fact be a general phenomenon in galaxies. 
Detections of other heavy elements in this absorber will require substantial improvements in the 
S /N of the spectrum (at the same resolution) . Perhaps the most promising transition to search for 
in future observations is the strongest Fe II transition at 2382.8 A. 

We have also estimated the metallicity of the Lya/O VI absorber at z a b s = 0.00337 in the 3C 
273 spectrum. We obtain [O/H] > —2.0 from photoionized as well as collisionally ionized models. 
If the gas is collisionally ionized, then the lower limit on the O VI/C IV ratio requires T > 10 53 K 
in the O VI phase. 

Several other sight lines to QSOs behind the Virgo supercluster have recently been observed 
or will be observed during the next year with the STIS E140M echelle mode. For example, it will 
be useful to examine the sight line to RXJ1230. 8+0115 (Read, Miller, & Hasinger 1998), which is 
only 54' from 3C 273. These nearby sight lines will likely provide additional valuable constraints 
on the heavy element enrichment of the IGM. 

We are indebted to Gary Ferland and collaborators for the construction of CLOUDY, and for 
making the code freely available. We also thank Jane Charlton, John Mulchaey, Bill Oegerle, Joop 
Schaye, and especially Ken Sembach for helpful comments. This research was supported by NASA 
Grant NAS5-26555. TMT acknowledges additional support from NASA LTSA Grant NAG5-11136. 
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A. Measurements of N(C II) and N(Si II) in High- Velocity Milky Way Gas in the 

Directions of HD116781A and HD122879 

As discussed in §4, there are relatively few Galactic interstellar clouds in which the Si II and C II 
column densities have both been reliably measured. This is primarily because N(C II) is difficult to 
measure in the Milky Way ISM; the C II resonance transitions are usually very badly saturated, and 
when the spin-changing C II] A2325 line is detectable, the corresponding Si II resonance transitions 
are strongly saturated. There is a spin-changing Si II] transition at 2335 A that can be observed 
when the C II] A2325 line is detected. However, detection of these semiforbidden transitions is 
challenging, and they have mainly been observed in sight lines such as ( Oph through relatively 
dense clouds. 

The recent high-resolution survey of interstellar C I absorption lines carried out by Jenkins 
& Tripp (2001) has fortuitously revealed a few additional high-velocity interstellar clouds in the 
spectra of two stars, HD116781A and HD122879, in which N(C II) and iV(Si II) can both be accu- 
rately measured. The continuum-normalized C II and Si II absorption profiles of these high-velocity 
features are shown in Figure 10. Due to the very high resolution of the spectra (i? ~ 200, 000), it 
is unlikely that these high-velocity lines are affected by unresolved saturation. Consequently, we 
can obtain reliable column density measurements. In the direction of HD116781A, at least two 
components are readily apparent in the high velocity gas. To deblend these components and obtain 
their individual column densities, we have used the profile- fitting software of Fitzpatrick & Spitzer 
(1997) with the line-spread function for this slit from the STIS Handbook. The heliocentric veloc- 
ities, Doppler parameters, and column densities obtained in this way are summarized in Table 4. 
Table 3 reports the [Si/C] values derived from these column densities. There may be evidence 
of additional components in these profiles, but these additional components are not sufficiently 
well-constrained to justify addition of more lines to the fits. We note that due to the weakness 
of the line, the wings of the Si II profile in the HD122879 spectrum might be lost in the noise. 
This could lead to underestimates of the b— value and N(Si II) in this cloud. Based on apparent 
column density ratios and fits to the Si II line forced to take on larger 6— values, we estimate that 
log iV(Si II) might be underestimated by 0.2 dex. The error bars in Table 3 reflect this source of 
uncertainty. 

A full discussion of the nature of these high- velocity features is beyond the scope of this paper. 
However, a few comments are in order. (1) The absorption lines occur at velocities well outside the 
velocity range expected for normal Galactic rotation in these directions (see Table 2 in Jenkins & 
Tripp 2001), i.e., these are "for bidden- velocity" clouds. (2) The [Si/C] abundances derived from the 
two high-velocity components toward HD116781A are quite similar (and consistent with the same 
value within the la uncertainties; see Table 3), but there are large differences in [Si/C] toward 
HD116781A and HD122879 at v Q « -109 km s" 1 (see Table 3). It seems unlikely that these 
differences are due to intrinsic changes in the interstellar abundances given the spatial separation 
between the sight lines (~190 pc at the estimated distance of HD116781A). As discussed in §4, it 
is possible that some of the differences in [Si/C] are driven by ionization. The differences between 
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the HD116781A and HD122879 sight lines are dramatic; 7V(C II) is only ~0.3 dex smaller toward 
HD 122879, but the Si II column density is reduced by an order of magnitude. Based on the modeling 
of Sembach et al. (2000), such large changes cannot be attributed to ionization effects if the gas 
is predominantly photoionized. However, this could occur if the gas is predominantly collisionally 
ionized (see Sutherland & Dopita 1993) and the gas temperature exceeds ~15,000 K. This would 
not be surprising since these clouds may well be shock-heated, given their large departures from 
normal Galactic rotation. In addition, some of the differences in [Si/C] between the sight lines 
may be due to more substantial dust grain destruction in the direction of HD116781A. This could 
occur if silicon is incorporated into more fragile portions of dust grains while C is locked into more 
resilient dust components. 
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Table 4. Forbidden- Velocity Absorption Lines in the Spectra of HD116781A and HD122879 



Sight Line Transition n (km s 1 ) a ' 6 6 (km s 1 ) a log iV a 



HD116781A.. 


C II A1334.53 


-107 


8.9 ±0.5 


13.68±0.01 




Si II A1260.42 


-109 


7.0 ±0.9 


12.10±0.02 




C II A1334.53 


-92 


5.9 ±0.5 


13.42±0.02 




Si II A1260.42 


-91 


3.6 ±0.9 


11.74±0.03 


HD122879.... 


C II A1334.53 


-109 


7.2 ±0.2 


13.40±0.01 




Si II A1260.42 


-110 


3.5l?; 5 7 


11.12±0.08 



a Doppler parameter, component velocity, and column density measured 
using the profile- fitting software of Fitzpatrick & Spitzer (1997). Uncer- 
tainties in this table only include random errors estimated by the profile- 
fitting code; systematic sources of uncertainty such as continuum place- 
ment are not included in these error bars. 

b For HD116781A, vlsr = v & — 5.9 km s _1 , assuming the standard def- 
inition of the Local Standard of Rest (Kerr & Lynden-Bell 1986). For 
HD122879, t> L SR = v© - 4.0 km s -1 . 
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(a) Milky Way sight line: HD116781A 



-150 -100 -50 

Heliocentric Velocity (km/s) 



Fig. 10. — Continuum-normalized absorption profiles of the C II A1334.5 and Si II A1260.4 lines 
in the Galactic sight lines towards (a) HD116781A, and (b) HD122879, plotted versus heliocentric 
velocity. The components of interest are indicated with tick marks. 
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Fig. 10. — (continued). 
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